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The authors introduce a geometry for ultrathin Au and Ag wires that ab initio calculations 
indicate to be more stable than previously considered planar geometries for these systems, by about 
0.1 eV per atom. This structure is insulating for both metals and for related Ago.sAuo.s alloys, with 
gaps of 1.3 eV for Au, 0.8 eV for Ag, and varying between 0.1 eV and 1.9 eV for the alloys. The 
insulating nature of the geometry is not a result of Peierls instabilities, and is analyzed in terms of 
an interplay between geometric and electronic structure effects. 
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Nanowires (NWs) based on 4d and 5d metals are a 
topic of intense current interest in the physics of nano- 
materials. Understanding the connection between atomic 
structure and transport, in the limit of the ultrathin 
monoatomic wires that have been produced experimen- 
tally, [l|, Q is crucial to the future manipulation of metal- 
lic wires and electric contacts in nanoscale electronic de- 
vices. Previous experimental and theoretical works on 
pure and alloy NWs have established that structure and 
transport are dynamically and strongly linked in the 
nanoscale limit. [2t3| Detailed investigation of the struc- 
tural landscape of ultrathin pure and alloy NWs is thus 
of paramount importance. 

Previous theoretical studies of infinite 

monoatomic Au and Ag NWs have considered two planar 
zigzag configurations, with angles between two noncolin- 
ear bonds along the chain of ^60° [see Fig. 1(b)] and 
^130°, respectively. Both zigzag geometries are stable 
for Au, while for Ag only the 60° geometry is a mini- 
mum of the E x i surface, pjjj This difference between 
Au and Ag has been tied to a stronger relativistic effect 
in Au, that leads to enhanced sd hybridization and stabi- 
lization of the low-coordinated structure, an effect that is 
also observed in the surface reconstruction of Ad and bd 
metals. [l3T | For both metals, the zigzag structures remain 
metallic, and their stability has been connected with a 
gapless Peierls transition (GPT), where transversal dis- 
tortions lower the energy, but a gap does not open at the 
Fermi level. In the recent work of Cheng and collab- 
orators, the pathway of the thinning process of Ag NWs, 
in mechanical break junction experiments (MBJE), 
is connected to the successive stable minima of the E x i 
curve, as the NW radius decreases. [l5| 

In the present work we examine the structure of ul- 
trathin noble-metal NWs. Our calculations indicate the 
existence of a planar structure for both Ag and Au, and 
also for Ago.5Auo.5 alloys, that has not been addressed 
in previous studies for these systems. In this geome- 
try, the average coordination remains fourfold, like in the 
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FIG. 1: (Color online) Ag, Au, and Ago.sAuo.5 nanowire 
structures: (a) SQ3; (b) ZZ4; (c) AiZZ4; (d) ZZ3+5; (e) 
AiZZ3+5; and (f) A 2 ZZ3+5. For the alloys, light (yellow 
online) and dark (blue online) circles represent Au and Ag 
atoms, respectively. 



60° zigzag one, but the four atoms in the unit cell alter- 
nate between threefold- and fivefold-coordinated sites, in 
a double zigzag structure [shown in Figs.[T][d)-(f)] that is 
lower in energy than the 60° zigzag. Moreover, the struc- 
ture is insulating, with energy gaps of 1.3 eV and 0.8 eV 
for Au and Ag, respectively, and does not result from ei- 
ther a Peierls or a GPT instability. Gaps for four possible 
related alloy geometries vary between 0.1 eV and 1.9 eV. 
No Peierls instabilities are expected for this structure, 
given its insulating nature. This is a striking result: the 
structure, besides being lower in energy than the planar 
ones considered in previous works, is an equilibrium (un- 
strained) insulating chain for these atoms, not related to 
Peierls instabilities, with sizeable energy gaps. 

We employ an ab initio methodology implemented 
in the SIESTA package, based on the Kohn-Sham 
formulation of density functional theory (DFT) . fl7T| 
The generalized gradient approximation (GGA) [18[ is 
used for the exchange-correlation energy, and norm- 
conserved Troullier-Martins pseudopotentials flfj| repre- 
sent the ionic cores. Cutoff constrained pseudo-atomic 
local orbitals form the basis of representation of the elec- 
tronic wave functions, with energy shifts of ~7 mRy. A 
double-zeta expansion is used for each radial basis func- 
tion, including polarization orbitals for the (I = 0) s func- 
tions. In all our calculations, convergence of total energy 
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TABLE I: Energy, length per atom, and electronic band 
gap (at the level of GGA-DFT) of Au, Ag, and Ago.sAuo.s 
nanowires, at the equilibrium geometry. 



1.0 1.1 1.2 1.3 1.4 1.5 
length/atom (A) 



1.0 1.1 1.2 1.3 1.4 1.5 
length/atom (A) 



FIG. 2: (Color online) Energy (in eV/atom) x strain (in 
A/atom) curves for Au and Ag nanowires in ZZ3+5, ZZ4, 
NPZZ4, and SQ3 geometries. 



differences and atomic forces to within 1 meV/atom and 
0.003 eV/A, respectively, was enforced. 

The structures of the infinite chains we consider are 
shown in Fig.[TJ (a) shows a planar threefold-coordinated 
square chain (SQ3). (b) shows the previously-considered 
60° planar zigzag fourfold-coordinated chain (ZZ4). (d) 
shows the planar double-zigzag chain we introduce in 
this study, with threefold- and fivefold-coordinated atoms 
(ZZ3+5). The latter is composed of two intercalated sub- 
units, as indicated in Fig. QJe) by primed and unprimed 
letters. Sites A and A' (B and B') are both threefold 
(fivefold) coordinated, and are related by an inversion 
symmetry in the ideal structure. We considered four alloy 
structures derived from ZZ3+5, with Au and Ag assigned 
to A and B sites, as follows, (i) AiZZ3 + 5 [Fig. EEfe)]: 
Au in A and A' and Ag in B and B'; (ii) A 2 ZZ3 + 5 
[Fig. EOT)]: Au in A and B and Ag in A' and B'; (iii) 
A 3 ZZ3 + 5: Au in A and B' and Ag in A' and B; and 
(iv) A 4 ZZ3 + 5: Ag in A and A' and Au in the B and 
B'. In addition, we considered two ZZ4 alloy structures, 
(i) A X ZZ4 [Fig. [He)]: with Au and Ag atoms alternat- 
ing on both the A (bottom) and B (top) sites; and (ii) 
A 2 ZZ4: with Au atoms on the A sites and Ag on the B 
sites. For all pure systems, we performed calculations at 
fixed strain (length per atom), with full relaxation of in- 
ternal degrees of freedom. For the minima of each E x I 
curve, full relaxation was performed, allowing cell vec- 
tors to relax until residual stresses were negligible. For 
the alloys, we only performed full relaxations to obtain 
the unstrained minima of the configurations described 
above. 

Figure [2] shows the E x i curves for each geometry, for 
Au and Ag. All energies are referred to the energy of the 
isolated atomic components. For Au (Ag), the minimum 
for the ZZ3+5 geometry is 0.12 (0.09) eV/atom lower 
in energy than the one for the ZZ4 structure. Among 
the structures in Fig. [2J the SQ3 geometry shows the 
highest equilibrium energy for both metals. Our results 
for the total energies, the equilibrium length, and the 
electronic band gap, for each structure, are included in 
Table I. For the ZZ4, our values are in good agreement 
with previously published results. @, f8l— fiol. \v\ 
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Regarding the stability of the ZZ3+5, we believe our 
results to be relevant in the context of the thinning pro- 
cess of the NWs in MB JE, as discussed in Ref. UH, when- 
ever the time scale of this process allows for the stabiliza- 
tion of the successive local minima, as the NW atomic 
density (and radius) decreases. In the atomic density 
range we investigate in this work, the ZZ3+5 geometry 
is the most stable, and should be a "magic structure" (in 
the language of Ref. Il5r i adopted by the NW in the thin- 
ning process. We can also look at the NW stability issue 
from the other end, i.e., at the ultrathin limit, by identify- 
ing possible transitions between locally stable NW struc- 
tures, as the applied tension, hence the atomic density, 
fluctuate. An upper bound for the barrier involved in 
the ZZ4-to-ZZ3+5 transition is given by the energy dif- 
ference between the ZZ4 minimum and the crossing of 
the ZZ4 and ZZ3+5 curves in Fig. [2j plus the barrier for 
the transformation, when both structures have the length 
per atom corresponding to the crossing point. In Fig. [2j 
we also show the E x £ curve for a non-planar variation 
of the ZZ4 geometry (NPZZ4), where small non-planar 
distortions are allowed. ZZ4 and NPZZ4 are very nearly 
degenerate at the corresponding minima. Our upper- 
bound barrier for the ZZ4-to-ZZ3+5 transformation is of 
the order of the room-temperature thermal energy for Ag 
(~50 meV) and Au (~70 meV). In the case of Au, a lower 
barrier is obtained starting from the NPZZ4, while in Ag 
a lower-barrier path starts from the ZZ4 itself. 

The band structures of the ZZ4 and ZZ3 + 5 structures 
for Au and Ag are shown in Fig. GHa)-(d). In both cases, 
the Fermi level (Ef) for the ZZ4 geometry does not occur 
at the it /2a vector, and the stability and semiconduct- 
ing nature of the period-doubling ZZ3 + 5 does not result 
from either a GPT or a Peierls instability. It derives, 
rather, from intricate band-structure effects, which are 
different between the two metals. We discuss first the 
band structure of the ZZ4 geometries. We have analyzed 
in detail the orbital-projected density of states resulting 
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FIG. 3: Energy bands for Au, Ag, and Ago.sAuo.s nanowires 
in ZZ4 and ZZ3+5 geometries. 

from the bands shown in Fig. |3] For both Ag and Au, 
the two bands crossing Ef are one almost-purely s band, 
crossing near one quarter of the Brillouin zone, and one 
hybridized sd band crossing near the point at three quar- 
ters of the zone, with predominantly d character. Due to 
the relativistic effect, the two metals differ in the nature 
of the hybridized band: s and d levels are more inter- 
mixed in the case of Au, while the s levels are higher in 
energy with respect to d levels, in the case of Ag. This 
explains the presence of a large pseudo-gap region below 
the Fermi level in the case of Ag which becomes very 
small in the case of Au, since the density of states (DOS) 
of the d-dominant band peaks at much lower energies 
with respect to the s band in the case of Ag. 

In the ZZ3 + 5 structure, the presence of two sites with 
different coordinations causes the DOS peaks dominated 
by the threefold (fivefold) sites to shift upwards (down- 
wards), relative to the ZZ4 peaks. In a tight-binding 
language, the onsite terms are shifted due the change in 
coordination. As a result, in the case of Ag, the gap is 
between two nearly purely s bands (plus the p character 
from the polarization orbitals). The top of the valence 
band (HOMO) is a bonding s band from the threefold 



atoms while the bottom of the conduction band (LUMO) 
is an antibonding band from the fivefold atoms. The 
bonding band from the fivefold sites is deeper in energy 
than the HOMO, and the antibonding one from the three- 
fold sites is higher than the LUMO. In the case of Au, the 
HOMO band is a strongly mixed sd band with d charac- 
ter from both types of sites (the fivefold site is dominant), 
and s character from the threefold sites, while the LUMO 
is the same as in Ag, the antibonding s band from the 
fivefold sites. 

For the alloys, the most stable AiZZ3 + 5 has Au in the 
threefold sites and Ag in the fivefold sites, while in the 
high-energy A4ZZ3 + 5 Au and Ag switch places. This 
is due to stronger relativistic effects in Au, that lead to 
its stronger tendency to form low-coordinated structures. 
The gaps for these two alloys can be obtained accurately 
from the HOMO and LUMO of the pure Au and Ag 
structures. Since HOMO bands are dominated by three- 
fold sites and LUMO ones by fivefold sites, and HOMO 
and LUMO energies are deeper in Au than in Ag, a larger 
gap is expected for AiZZ3 + 5, with Au at threefold sites 
leading to a deeper HOMO, and Ag at fivefold sites lead- 
ing to a higher LUMO. The opposite effect would occur 
in A4ZZ3 + 5. This is indeed what we obtain, as shown 
in Table I. Further, Table I shows a clear correlation be- 
tween energy and band gap for the alloy structures. Note 
that the insulating AiZZ4 alloy has a lower energy than 
the A 4 ZZ3 + 5. 

In summary, the authors introduce a geometry for ul- 
trathin Au, Ag, and Au-Ag-alloy wires that is more stable 
than previously considered planar geometries for these 
systems, by about 0.1 eV per atom. The insulating na- 
ture of this structure, with gaps of 1.3 eV for Au and 
0.8 eV for Ag, is not related to Peierls instabilities, re- 
sulting, rather, from an interplay between geometric and 
electronic structure effects. 
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